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A PREDICTION MODEL FOR LIFT-FAN SIMULATOR PERFORMANCE 


By 

Joseph A. Yuska 

Advisor: Dr. Chittaranjan Jain 


ABSTRACT 

In' this report the' performance characteristics of a mode] VTOL 
lift-fan simulator installed in a two-dimensional .wing are pre- 
sented^ The lift-fan simulator consisted of a 15-inch diameter fan 
driven by a' turbine contained in the fan hub. The performance of 
the lift-fain ‘simulator was measured in two ways:' (1) the calculated 
momentum thrust of the fan and turbine (total thrust loading), and 
(2) the axial^-force measured on a load cell force balance (axial- 
force loading). ' Tests were conducted over a wide range of crossflow 
velocities., corrected tip speeds , and wing angle of attack. 

The objective of this thesis work was- to develop a prediction 
modeling technique to help in analyzing the performance characteris- 
tics of lift-fan simulators. A multiple linear regression analysis 
technique is presented which calculates prediction model equations 
for the dependent variables: total thrust loading and the axial- 

force loading as a function of the independent variables: crossflow 

velocity, corrected tip speed, and wing angle of attack. The re- 
sults of the regression analysis show that there was a significant 
correlation for the two loading parameters as a function of the 

ii 



independent variables. The results also showed that there is a sig- 
nificant correlation for the total thrust loading as a function of 
the axial-force loading. 

The test results showed that the effect of increasing crossflow 
velocity is not the same for total thrust loading as for axial-force 
loading. The multiple linear regression analysis was used to deter- 
mine a significant correlation for the difference between the, two 
loading parameters for the same independent variables.' Crossflow 
velocity was found to be the most significant variable in the force- 
loading difference correlation equation. 

The results of this 'work indicate that, this prediction modeling 
technique could be valuable in designing and -analyzing lift-fan sim- 
ulator tests. It could be used to reduce the amount of instrumenta- 
tion required and' tp reduce the number of test runs required to de- 
fine the performance characteristics- of lift-fan simulators. The 
Use of this technique could reduce the cost of building and testing 
lift-fhn simulators. 
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1.0 INTRODUCTION 


Vertical take-off and landing (VTOL) aircraft are currently 
under study as a means for improving short-haul intercity air trans- 
portation systems. VTOL aircraft Can relieve airport congestion,, 
reduce air time delays, and can also supply access to communities 
currently without air transportation. The feasibility of VTOL in- 
tercity air transportation systems has been explored by a number of 

‘ ’ j. 4* 

domestic and foreign organizations and government agencies ’ • It 

• is recognized that there are many aspects to the problem of VTOL 

transportation system - economic, sociological, and political, as 

well as technological^.; 

The' bulk of the -research oh VTOL transportation systems has 
been conducted on the technological problems. Some of the. earliest 
. wind tunnel test of lift fan concepts were performed in 1957 by 
Hickey for the National Advisory Committee for Aeronautics .of the 
United States^** 7 , Wardlaw, Templin and McEacherii for the National 

Aeronautical Establishment' of the. National Research Council of' 

8 9 • 

Canada ’ ,'and Gregory and Raymer for .the National Physical Labor- 

• ' • 10 
atory for the Ministry of Aviation, London ,• England . These tests 

were performed on fans and propellers having low pressure ratios. 

They found that the lift and drag forces and pitching moment on the 

model could not be predicted by considering only the force reactions 

due to the fan momentum thrust and the wing lift. It was surmised 

* 

' Numbers refer to the references listed at the end of the report. 
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that part of these discrepancies could be accounted for by consider- 
ing the changes in the flow field about the wing due to fan-wing 
interaction. 

All three countries have continued research to determine the 
fan-wing-flow interactions. In England, the investigations were 

carried on with scale models having low fan-to-wing area ratios. 

/. 

Results of this work have been published in references 11 through 

.14 

■ 13; ■ In Canada, tests were also performed on scale models . . In the 

United States the bulk of the work was carried out at Ames Research 

*15 16 

Center with large scale models 5 . Hickey (in reference 16) has 

proposed a rudimentary method for estimating the lift, drag and 
pitching moment which gives good agreement with his ’ experimental • 

work. Other NASA work on fan-in-wing concepts was performed by 

• 17 18 
Mort and Spreeman 

19 • • .. • . 

Norland attempted to - develop a . two-dimensional theory to ex- 
plain the decrease in lift- at low forward speeds which has-been in- 
dicated, in small scale tests of fan-in-wing configurations. -Correl- 
ation of Norland's theory with his experimental work is inconclu- 
sive, however it provided a basis for a sound theoretical treatment 
of the fan-in-wing problem. 

• Feasibility studies of V/STOL concepts for 'short-haul transr . 
•pprts has .shown the need for advanced-concept lift fans that are 
quiet, smokeless, and reliable. To develop these advanced-concept 
lift fans requires detailed ‘design data on the performance of 
advanced-concept lift fans in the transition speed range (crossflow 
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velocity) from vertical to normal wing-supported flight, There are 
a number of research problems associated with the performance of 
lift fans operating in the transition speed range. One problem, 
crossflow inlet distortion encountered during transition, has an 
adverse affect on lift-fan performance and noise generation. A 
second problem is the effect of ' crossflow velocity changing the 
back pressure at the exit of the lift fan which changes the opera- 
ting point of the lift fan. A third problem is the interference 
effect of the interaction of the fan flow and external airflow. 

NASA Lewis Research Center recently embarked on a research 
program to generate detailed design and performance .data on ad- 
vanced lift-fan concepts.' An initial experimental test project 
measured the detailed crossflow performance of a 15-in. diameter 
lift-fan simulator installed in a two-dimensional wing. The fan- 

in-wing model tests were conducted in a 9- by' 15-ft V/STOL wind 

* ,20 • ■ • ' • ' 

tunnel .• ... 

Results of the measured fan axial force- data and overall aero.- 

dynamic characteristics of the fan-in-wing configuration have been 

‘21 ' • 

reported by Yuska and Diedrich . The -applicability of these re- 
sults to the thesis topic will be discussed in other sections of 
the thesis . 

Presently a report is in preparation to present the internal 
flow characteristic of the lift-fan simulator, and analysis is 
being conducted to determine the relationship between the axial 
force measured on the load celi force balance and the calculated 
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momentum thrust of the fan and turbine. In addition, tests have 
been completed on three other lift-fan simulator designs. 

Since lift-fan simulator models are very expensive and test 
time is limited and expensive, it is desirable to develop a mathe- 
matical prediction model for the lift-fan simulator performance 
characteristics over its entire operating range. 

The purpose of this thesis is to present the results of a mul- . 
tiple linear regression analysis prediction model • for a 15-in.- 
lift-fan simulator. The prediction model is developed from ob- 
served values. of lift-fan simulator thrust loading and axial-force 
loading as the dependent variable and major test variables (cross- 
flow velocity, fan corrected tip speed and wing angle of attack) as 
the independent variables. 

There are a number of ways that the prediction model may be 
used. A reliable prediction model' could be used to reduce the num- 
ber of test runs required to define the performance, characteristics 
of the lift-fan simulator. 'The use of the prediction- model would 
allow the experimental investigator to make better use of the avail- 
able test time and reduce the cost of experimental test programs. 

The use of replicated data could provide an' estimate of the accuracy 
to the test data. The prediction model could also be used in the 
design of future experimental lift fan simulators. • 

The prediction model would also be useful in smoothing the data 
for the analytical work of determining the relationship between the 
axial force measured, on the load cell balance and the calculated ' 



- 5 - 


momentum thrust of the fan and turbine. This relationship is use- 
ful in the analysis to determine the interference effect of the fan 
flow and external airflow on the lift-fan simulator. 



2.0 DESCRIPTION OF EXPERIMENT 

The objective of the experiment was to determine the perform- 
ance of a model VTOL lift-fan simulator installed in a two- 
dimensional wing section for various test conditions. The lift-fan 
simulator consists of a 15-in. diameter fan driven by a turbine con- 
tained in the fan hub. The wing section was designed to simulate 
the aerodynamic flow which 'would be produced by a VTOL transport 
aircraft during landing: or take-off. 

• The performance of the lift-fan simulator is measured by the 
total axial thrust developed under various conditions. The. total 
axial force, F, is the sum of the fan momentum thrust, Fp, lb; the- 
.turbine momentum -thrust, F^, lb; and the fan base forcej Fg, lb or. . 

F = Fp + Fp + Fg (l) 

The fan and turbine momentum thrust is calculated from observed 
values of total pressure, P, total temperature, T, and. airflow’ • 
angle, .(3, ■ at. the fan and turbine exit; and from observed 'values of 
static pressure, p^, at the fan inlet which is correlated with fan. 
weight flow, Wg, using a correlation equation from a theoretical 
technique developed by Stockman^, and from observed values of tur- 
bine weight flow, wm, measured with an orifice plate. Another 
measure of fan performance was the axial force F^q measured on a 
load cell force balance., .The 'force balance consists of three 
500 lb load cells equally spaced on a 22.5-in. diameter circle.. 


-' 6 - 
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The load cell force balance attaches the lift fan simulator to the 
wing. , x : 

The major test variables were: 

(1) Free-stream (airflow) velocity, V, feet per second (this, 
parameter is also commonly called crossflow velocity in VTOL appli- 
cations and may be used interchangeably). 

(2) Fan corrected tip speed, U^/t/e" expressed as a- percentage 
of the corrected tip speed at design conditions (U^/v^T = 980 ft/sec 
(Tip speed is the product of the angular velocity of the fan rotor 
blade times the tip radius of the rotor blade; and 0 is the ratio 

a ’ * . 

of inlet total temperature to standard sea-level temperature of 
518.7° R). 

(3) Wing angle of attack, a, deg. 

Tests were conducted at many possible combinations of test var- 
iables. However, the test conditions can be grouped into two major 
groups:- (1) Static test • conditions where the crossflow velocity is 
zero, and (2)- crossflow test conditions where the crossflow ' velocity 
varied from’ 70 to 240- feet per second. For- static test conditions, 
test data were observed at eight values of fan corrected tip speed 
and two angles of attack. - For crossflow test conditions, test data, 
were observed at six crossflow velocities, two values of fan correc- 
ted tip speed, and eleven values of angle of attack. Test data were 
observed at 52 points . corresponding to different combinations of ' 
test points. Table I is a summary of the observed test conditions. 
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TABLE I 

SUMMARY OF TEST CONDITIONS AND INPUT DATA 


Crossflow 

Corrected 

Angle of 

Total 

Axial- 

Force 

velocity. 

tip speed, 

attack. 

thrust 

force 

loading 

v, 

Uij* / y^j", 

a > 

loading , 

loading , 

difference 

ft/sec 

• percent of 

deg 

F/<5A, 

F/6A, 

&F/6A, 


design 

.speed. 


lb/ft 2 .. 

lb/ft 2 

lb/ft^ 

82.1 . 

100; 38 • 

. -10.10 . 

758.2' 

• 701.7, 

-56.55 

117.5 

99.91 

-10.11 

731.4 

- 691.2 

-40.15 

151.0 

99.20 

‘ -10.11 

709.9 

660.3 

-49.52 

180.0 

99.88 

-10.29 . 

702.6 

667. 0 

. -19. '90 

.207.4 

99.87 

-10.12 ■ 

■ 692 ; 0 

672.0. 

-35.58 

239; 9 

99,61 

-10.16 

6-90.4 

. 677.1 

-13.30 

76.8 

100,18 

-0.22 

. 760.9 

716.3 

-44.55 

. 120.1 

99.74 . 

-0.20 

731.0 . 

. 709.8 

-21.17 

151.9 . 

.. ,95.47 • • 

-0,17 

... 711.3 

.682.4 

-28.92 

180.2- 

99 ! 34 

-0.24 

698.9 

687.9 

-11.03 

203.3 

99.90 

-0.15 

700.1 

722.0 

21.95 

236.4 

99.21 

-0.20 

679.2 

726.8 

47.57 

71.8 

100.06 

10.06 

764.0 

721.1 

-42.85 

. .118.3 

99 ,57 • 

’ 10.06 

739.5 

■725.6 

-13.93 

153.0 . 

99. .53 ' 

• 10.03 ' 

715 .3 

'728,2 

12.95- 

177.8, 

,99.54 

. 10.04 

699,2 

734.6 

•35.44 '. 

,204.1 

99.90 . 

10.03 . 

680.6- 

' 747.1 

.. 66.49 

234.4 

• 100.66 

10.10 

'649.9 

756.6 

106.73 

81.3 

70.34 

-10.10 • 

376.6 

' 343.7 

. -32’. 86 

118.0 

‘ 69.78 

-10.10 

375.8 

333.6 

-42.23 ’ 

158.6' 

70.01 . 

-10.20 

378.8 

342.6 

-36.20 

181.8 

. 69.74 • 

. -10; 20 

380.3 

. 357.4 ' . 

-22.92 ■ 

,205.1 

69.69 

. -io.io 

377.8 . 

370.3 

-7.53. 

,' 238.7 

68,88 

-10.20 

. 382.0 ' 

391; 2 

9.21 

78,0 

70.27 

-0.20 

380.3 

•346.2 • 

■ -34.07 

121.8 

70.25 

-0.20 

378.7 

357.6 • 

-21.08 
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TABLE I (Concluded) 

SUMMARY OF TEST CONDITIONS AND INPUT DATA 


Crossflow 

Corrected 

Angle of 

velocity, 

tip speed, 

attack , 

V, 

Uj/ , 

a. 

ft/sec 

percent of 
design 
speed 

deg 

151 . 3 

69'. 78 

-0.30' 

178.5 

69.85 

-0.10 

206.4 

69.40 

-0.20 . 

237.2 ' 

69.23 

-0.10 

74.0 

'70.41. • 

10.10 

120.5 

70.20. 

10.10 

153.3 ' 

'70.06 • 

' 10.20 

" 182.1 

69.61 

• 10.20 

204.5 

69.82', 

10.20 

234.3 

. 69.70 

•10.10 ' 

205.6 

69.46 

-13.60 

204.4 

69.39 

-5.10 

207.5 

69.73 

4.90 

204.4 

69.84 • 

15.50 

204.2' 

99.75 

• rl5. 30 

204.4 

99.93 

-5.10 

201.5 

100-, 23 

5.00 

204.6 

99.15 

15,50 

0 

70.83 

0 

0 

• 81.14 

0 . • 

0 

85.89. 

0 . 

0 

91.57 

0 

. 0. 

.96.28 

. 0 

• 0 

101.39 

0 

0 

106.18 

0 

0 

111.73 

0 . 


Total 

Axial- 

Force 

thrust 

force 

loading 

loading , 

loading, 

difference 

F/6A, 

O 

F/6A, 

O 

AF/<SA, 

O 


lb/ft 2 lb/ft 2 ' lb/ft 2 


381.3 

372.8 

■ -8.50 

379.1 

387.3 

8.19 

385.5 

412.3 

26.78 

375.0 

417.3 

42.35 

381.3 

365.0' • 

-16.33 

378.4 

376.8 

. -1.55 

380.1 

395.6 

•' 15.46 

381.2 

417.0' 

. 35.82 

371.5 

437.6 

• 66:14 

345.0 

451.9 

106.. 86. 

381 ! 4 

355.9 

-25.53 

379.7 

386.8 

7.06 

377.1 

' 434.1 

• 56.99 

348.3 

440.6 

92.30 

692.3 

643.9 

' -48.20. 

694 . 5 

686. 0. 

-8.. 50 

691.5 - • 

726.9 

•35.40 

66-7.8 

■ 766.0 

. 98.20 

411.4 

. 347.9 

-63.50 

538.7' 

• 479.8 '■ 

-58.90 

613.3 

'542.4 

-70.90 

696 ,0 

624.7. 

• -71.30 

762.6 

691.7 

.-70.90 

819.2 

745.8 

-73.40 

860.1 

•802 . 3 

• -57.80 

910.1 

852.9 

-57.20 



3.0 EXPERIMENTAL MODEL AND APPARATUS 


In this chapter, a' detailed description of the experimental 
model and apparatus used in this experiment is presented. The test 
procedure, instrumentation, and data recording equipment used for' 
data collection are described. 

3.1 Wing and Wing Balance 

A photograph of the fan-in-wing model installed in the 9- by 
15-ft V/STOL test section is shown in Figure- 1. The wing spanned 
the 9-ft height of the test section and had a 4.5-ft chord. The 
wing has a constant NACA 65^(218) -217 airfoil section along, the 
span except for localized modifications required' for fan installa- 
tion. 

The fan-in-wing model was attached to a wing balance that 
measures overall model lift,' drag, and pitching moment. Figure 2 
is a sketch of the. wing balance showing pertinent dimensions.' 

3.2 Fan and Drive Turbine Characteristics 

The fan assembly was; attached to a .force balance System within 
the wing. The fan force balance consists of three 500-lb load cells 
equally 'spaced - on a 22.5-in. diameter circle, as shown in Figure 3..- 
The balance measured- the combined axial force- of the fan and turr • ' 
bine. Drive air was supplied from both sides of the wing. 

Figure 4 shows a cross section of. the model lift fan. The fan 
rotor tip diameter was 15.2 in., and the stage exit diameter was 
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Figure 1. - Fan-in'-wing model in 2. 74- by 4.58-meter (9- by 15-ft) V/STOL . 
wind tunnel. . 
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14.78 in. The diameter of the fan inlet bellmouth at the wing 

junction was 22.0 in. The fan inlet bellmouth was contoured flush 

with the wing upper surface. The fan hub protruded from the wing 

lower surface as shown in Figure 4. The inlet and exit flow area 

2 2 

of the fan stage were 0.990 ft and 0.868 ft , respectively. The 
fan rotor had 35 blades which were made from 18 Nickel-200 maraging 
steel. Two 0.060-in, diameter damper wires located at 60 percent 
. of the passage height from the hub provided sufficient damping to 
prevent blade flutter and reduce blade vibration. . • 

The fan drive was a compact, two-stage supersonic turbine con- 
tained in the fan hub. The turbine was scaled from a previous NASA 

23 24 * 

design * -. . .The turbine had a mean diameter of .5.4 in. and devel- 

oped 685 hp to drive the fan at design conditions. The turbine was 
driven by high-pressure air supplied through- six equally-spaced 
rstruts (1'2.5 percent thickness ratio) spanning the fan -passage 
(Fig. 4)-. The struts connected - the inner and oater 360° plenums . 

For some tests, exit louvers were mounted, on the wing as shown 

in Figure 5. For this installation, ‘the louver forces appeared 

\ * 

only on the wing- balance and not on the fan balance- • .'The exit lou- 
vers covered -both the fan and drive turbine exhaust. The exit lou- 
vers were capable of remotely vectoring. the fan and drive turbine 
thrust from 3° forward (negative louver angle of attack) to 40° aft- 
(positive louver angle of attack) . . 
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Figure 5. - Exit louver installation. 
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3.3 Instrumentation for Data Collection 

The major features .of fan internal instrumentation are shown 
in Figure 4. Static pressure taps were located at eight circumfer- 
ential positions on the outer shroud of the fan inlet bellmouth 3/4 
of an inch above the fan rotor leading edge. These eight static 
pressure measurements were averaged together to provide an average 
inlet static pressure for use in the weight flow correlation.- At 
the rotor exit, six equally spaced stator blades were each instru-. 
mented with total temperature -probe's at three radial positions. 

At the duct exit, six fixed-position rakes equally spaced' 
around the circumference- (positioned between the turbine drive-air 
struts) -were used to measure, the airflow angle at three radial po- 
sitions and total pressure at five radial positions as shown in 
Figure 6. Static pressure taps were located at six equally-spaced 
positions on the -hub an'd tip Surfaces in the measuring plane of the 
total pressure probes. The airflow angle- probes were fixed- 
position, non -nulling -instruments . The construction details and. " • 
calibration techniques for fixed-position flow direction probes are 
reported' in reference 25. ' 

The turbine exit duct was instrumented at two circumferential 
locations to measure flow angle .at one radial position,, total pres- 
sure at five radial positions, and total temperature at two radial 
positions (Fig. 7). Two equally-spaced static taps were located on 
the turbine exit inner and outer wall surfaces. Turbine air weight 
flow was measured by a calibrated flat plate orifice. Base static 




Figure 5.. - Fan Duct Exit Plane Instrumentation Location. 


Hub static 



Figure 7. •- Turbine Exit Plane Instrumentation Location. 
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pressures were measured on the ring between the fan and turbine 
flow passages and on the centerbody surface of the turbine. These 
static pressures were used to calculate the fan base force. 

A summary of the total instrumentation provided for in the 
experiment .is shown in Table II and Table III. A total of 313 
pressures and 31 temperatures were measured. To reduce the number 
of strain gage, pressure transducers required to'. measure the 313 
pressures, a scanning type pressure sampling valve for measuring . 
multiple pressures was used. ' . 

Each valve is recorded in- sequence and controlled by. the data 
recording system. • The valve rotor is then stepped to' the next port 
for the next measurement.- The process is repeated until all pres- 
sures are measured. The stepping rate may be varied to match the 
settling out time for the pneumatic pressure system. The trans- 
ducer signal is connected to the automatic voltage digitizer and to 
the data recording system. With this systemone pressure trans- • 
ducer can measure 46 individual pressure's. In this experiment 
seven pressure sampling values were used. Each .using only one pres- 
sure transducer. Thirty two additional transducers were used to 
measure various required gage and absolute pressures, and pressures 
beyond the measurement range of the transducers in the sampling val- 
ves. Thus, only 39 pressure transducers were used to measure 313 
pressures. ' . 



TABLE II 


SUMMARY OF TOTAL NUMBER AND TYPE OF PRESSURE 
AND TEMPERATURE MEASUREMENTS 


1 

Location of instrumentation 
in lift-fan simulator 

Total number 
pressure 
measurements 

of 

Total number of 
temperature 
measurements 


P 

P 

AP 

P A 

T 

AT • 


s 



A 



Wing surface 

24 •' 




3. 


Bellmduth. front 

40 






■ Bellmouth. rea.r .. 

24 

12' 





Rotor exit plane 

12 

54. 

18 

18 '! 


•• 18 . 

Duct exit plane 

12- 

• 30 

18 

18 



Turbine inlet orifice ■ 

1 


3 


• 1 


■ Turbine inlet header 

1 




3 


Turbine . inlet plenum 

' 1 






Turbine exit plane 

• 4 

. 10 

2 ' 

2 

. .. 4 


Fan .exit lower Surface 

3- 





. 

Tunnel conditions 

3 


3 


'2 


TOTAL 

' 

125 ’ 

106 

44 

'38 

•13 

18 


' . ' TABLE '.III .. . 

. SUMMARY OF TOTAL NUMBER AND TYPE -OF -FORCE AND 
MISCELLANEOUS DATA MEASUREMENTS 


..Type of measurement 

Total number of 


' measurements 

Lift force . 

3 

Drag force 

2 

Axial force 

3 

Fan RPM- 

1 

Fan torque 

1 .. 

Wing angle of attack 

1 

Louver position 

1 

Fan box- static pressure 

1 

Rotor strain 

2 






3.4 Data Recording Equipment 

The data is recorded on the Lewis Research Center central data 
recording system called the Central Automatic Digital Data Encoder 
(CADDE) that is used by many test facilities to record digital read- 
ings' from transducers of pressure, voltage, and events per unit 
time. A schematic of the system is shown in Figure 8, and reference. 
26 gives a detailed description of the system. The system has a 
300-rdata. Word capacity. In this experiment 368-data words were re-- 
corded at. each test point. The data are recorded on magnetic tape 

as binary-coded decimal digits with additional characters for idep- 
. * . . . 
tification and computation instructions. The magnetic tape is the- 

'permanent record of the raw data. In .addition, the data may be . 

simultaneously recorded on an IBM 360 time-sharing computer system. 

This allows the test engineer to receive "on-line" computer data • 

typed out in the control room after the data point, is taken. After 

recording, d duplicate tape- is made .to process the data as a - stored 

program on' an IBM 709-4 computer. 

3.5 Experimental Test Procedure 

The general test procedure- was to set the desired. test condi- 
tions of crossflow velocity, corrected tip speed; wing angle of 
attack, and exit louver deflection angle. The fan internal • flow 
and tunnel flow conditions were allowed to stabilize at a steady 
state condition before the start of .each data recording point. The 
time to acquire data (data scan' time) for a single test point was 
typically about 25 seconds.. During the data acquisition period, 
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test conditions were stabilized through the individual control sys- 
tems for each of the preceding independent parameters. The fan 
speed generally fluctuated somewhat during the data acquisition 
time. Over the speed range, the variation was +50 rpm which was 
only a +0.3 percent variation of the design speed. 

Test Procedure for Static Test Conditions (Zero Crossflow Velocity) 
The tunnel circuit was vented to the atmosphere to' minimize 
recirculation of the fan exhaust. Static tests were conducted with 
the fan axis both perpendicular and parallel to the tunnel axis. 

The orientation with the fan axis parallel to the tunnel .axis and 
the inlet of the fan facing upstream was used to achieve the most 
steady flow conditions at the fan inlet. Data taken with the fan 
axis perpendicular to the tunnel axis showed that the internal flow 
of the fan was less steady than with the fan axis oriented parallel 
to the tunnel axis, although the: overall thrust levels were' the 
same for both orientations. * 

Test, Procedure for Crossflow Test Conditions 

The general test procedure was to first bring the tunnel up .to 
speed and set the desired level of crossflow velocity in the test 
section. Before the start of each data recording point,, the tunnel 
flow was allowed to stabilize for several minutes while the desired 
values of wing angle of attack, louver position, and fan speed were 


established. 




4.0 DATA ANALYSIS TECHNIQUES 


This chapter describes the data analysis procedure used in de- 
veloping a prediction model for the lift-fan simulator. The de- 
pendent variables, total axial force, F, and the axial force meas- 
ured on the load cell force balance, F^, were calculated for the 
various test variables. The details of these calculations are pre- 
sented in Appendix A. To enable the test results of various lift 
fans to be compared, the total axial 'force and the axial force 

measured on the load cell- balance are- divided by the fan frontal 

2 

area, A = 1.263 ft ', and 6 the ratio of inlet total pressure to 

. * • 2 

standard sea- level pressure of 2116.2 lb/ ft . The resulting ratio, 
F/SA or F tp /< 5A is the total thrust loading parameter or the' 
axi^l-force loading parameter, respectively. These parameters nor- 
malize the data so that the results from this lift-fan simulator 
test can be compared to the test results of other lift-fan simula- 
tors or full-scale lift fans having different fan' frontal areas. 

The calculated dependent variables and the independent test 
variables were then used as input data (shown in Table I) to deter- 
mine the lift-fan simulator prediction model. A multiple linear 

27 

regression computer program (RAPIER) , which is a program in the 
NASA Lewis Research Center computer program library, was used to 
compute the lift-fan simulator prediction model equation. A de- 
scription of the multiple linear regression computer program is 
presented in Appendix B. 


25 
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The model equation was assumed to be 

y = b Q + b x v + b 2 (u T //e) + b 3 a + b 4 v(u T //e) 

+ b 5 Va + b 6 (U T //0)a + b 7 V(U T //0)a 
+ b g V 2 + b g V 2 a + e (2) 

where y is the dependent variable corresponding to F/6A or 

F t „/ 6A; b«, b ,*<*b D are. the regression coefficients and e is 

lo u i. y 

the error term. The computer program was then used with the back- 
ward rejection option to delete regression coefficients that have 
t-statistics which. are less than 90 percent significant. The pro- 
gram then computes the final regression coefficients, the analysis 
of varience (ANOVA) table and the coefficient of correlation, r. 

An F test is made for testing the null hypothesis H^: that 

there is no correlation, that all the correlation coefficients are 
zero. The . alternative hypothesis is H^: that there is a significant 
correlation. If the calculated value of F is greater than the 
critical value of the F distributions for the. number of degrees 
of freedom and the specified significance levels, then the null hy- 
pothesis Hq is invalid indicating that there is a significant 
correlation and the regression coefficients are not zero. The pro- 
gram then computes a calculated response (y) for each set of inde- 
pendent variable observations using the prediction model equation. 

A comparison is then made between the observed values and the cal- 


culated values. 
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Since there was an observed difference between the total axial 

force loading and the axial-force loading measured on the load cell 

balance, another dependent variable was calculated for all test 

variables (see Table I). This variable is called the force loading 

difference, AF/6A, which is the difference between F /6A and 

LC 

F/6A. 


AF/6A = F T _/6A - F-/6A (3) 

The multiple linear regression computer program was again utilized 
to determine if a correlation existed for force loading difference . 
as a function of the independent test variables. " . . 

The RAPIER computer program also 'has a prediction option, 
which used the regression equation to predict an estimated response... 
at some condition of the independent variables. This option is use- 
ful in calculating responses at fixed values of the independent var- 
iables. This technique is useful in making carpet plots to allow 
the presentation of the dependent variable as a function of .all its 
independent variables on a single graph. 



5.0 RESULTS AND DISCUSSION 


la this chapter, the results of the multiple linear regression 

analysis are presented. The correlations were obtained for total 

thrust loading, axial-force loading, and the difference between the 

axial-force loading and total thrust loading. In each of the above 

cases the independent variables were crossflow velocity, corrected 

tip speed, and angle of attack. In addition a final correlation 

was obtained for the total thrust loading as a function of axial- 

force loading. The results of the statistical computations are. 

discussed, and a comparison of the correlation equations and the 

% • 

observed vaiues is presented. ■ The physical significance of the 
observed variations in loading parameters are also interpreted for 
the model VTOL lift-fan simulator. 

5.1 Correlation • Equations 

The prediction model correlation equation for each dependent 
variable was computed using the input data shown in Table I. The 
resulting correlation equations are shown in Table IV. From Table 
IV it can be seen that in each correlation equation for total 
thrust loading ,. axial-force loading, and force loading difference, 
certain regression coefficients were deleted by the backward rejec- 
tion option of the computer program because the t-statistics test 
of each delected regression coefficient showed a significance level 
of less than 90 percent. 
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TABLE IV 

CORRELATION EQUATIONS 


1. Total thrust loading: 

F/6A = - 529.2807 + .8361 V + 13.1970 (U rj? / /e") 

- .0132 V(U T //0) + .01586 Va - .00010 V 2 a (4) 

2'. Axial-force loading: : • . 

F / <5A = - 566.0920 + ..9242 V + 12.9196 '(.U_//e) 

LG i. 

• ‘ . - .0132 V(U T //0) + .02823 Va- - .02678 (U /*^0)a • 

+ .000268 V(U T //0)a + .00158 V 2 - .00011 .V 2 a (5) 

3. Force loading difference: 

AF/6A = - 38.7907 + .1294 V - .2679 (U T //e) 

. . + . 00007 9 V (U^/Zq ) ql + .00141 V 2 

+ .00007 V 2 a ' • . • ' • (6) 

4. Total thrust loading as a function of axial-force loading: 

F/6A = 2.697+ 1.016 F / 6A 


.(7) 
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Equations 4 and 5 in Table IV show that velocity V and tip 
speed U^,/ are the dominant independent variables in the equa- 
tions for total thrust loading and axial-force loading. In both 
equations the regression coefficients associated with V and 
U,j,//0 are positive and have nearly the same value. The equations 
for total thrust loading (equation 4) also contains three interac- 
tion terms (V(U '//0), Va and V a). The equation for axial-force 
loading (equation' 5) contains six interaction terms (V(U^,//o"), Va., 
(U^,//6)a, V(U T //&)a, V , and V a). These interaction terms ac- 
count for the observed variations (increases and decreases) in 
total force loading and axial-force • loading which will be discus- ■ 
sed in a later section of this report. 

Equation 6 in Table IV shows that velocity V and tip speed 

are again the dominant independent variables in the force loading 

difference equation. The equation also contains three interaction 
, — 2 2 

terms .(V(U^/v0), V , and V a). The regression coefficients of the 
interaction terms are all positive which indicates that at a con- 
stant tip speed as V and a increase the value of AF/6A will 
increase. The observed increase in AF/6A with increasing V. and 
a will be discussed later. The correlation .equation for total 
force loading as a linear function of axial-force loading is shown 
in equation 7 of Table IV. ’ ■ 

For each regression equation, the correlation coefficient was 
calculated: and the null hypothesis H^: that there is no correla- 

tion, was tested. An F-test was performed, where the calculated 



F value is compared to the critical value of the F distribution 

obtained from Table 7.9 in reference 28 for a 0.01 probability of a 

larger value of F and the computed number of degrees of freedom. 

The critical value of F is denoted as F_ n . where the 

0 ’ 01 > v 2 » V 1 

subscript 0.01 is the specified probability level; and \ is the 
degrees of freedom for the mean squares of the regression; and 
is the degrees of freedom .for the mean squares of the residual. 

The results of the Analysis of Variance are shown in Table V. In 
each correlation equation for , total, thrust loading, axial-force 
loading, and force loading difference, the F-test showed that for 
each of these three correlation equations, the computed value of F 
was greater than th'e critical value of the F-distributions for the 
specified probability level of 0,01. Therefore, the null hypoth- 
esis HqI that there is no correlation, that all the coefficients 
are zero, was rejected; and the alternative hypotheses H^: that 

there is a significant correlation* was accepted.. In each case, 
the correlation coefficient was greater than 0.988 which again in- 
dicates that there is a significant correlation for each dependent 
variable. 

The correlation equation for total force loading as a function 
of axial-force loading was assumed to be -linear. This assumption 
would probably be valid for zero crossflow conditions but not for 
nonzero crossflow conditions. It has already. been seen that the 
effect of crossflow velocity on both loading parameters is non- 
linear. Previous experience has not shown how the total thrust 
loading varies as a function of axial-force loading. Even with the 
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linear assumption, the correlation equation for total thrust load- 
ing as a function of axial-force loading has a high degree of cor- 
relation. The correlation coefficient r is 0.9620, and the over- 
all F-test indicated that the alternative hypothesis H^: that 

there is a significant correlation, should be accepted. 

5.2 Comparison of Correlation Equations and Observed Values' 

The values calculated from, the correlation equation and the 
observed values of the dependent variable for the various test var- 
iables are compared. The comparison of the calculated and observed 
values for the total thrust loading and the axial-force loading are 
shown in Figures 9 thru 12. The data on Figure 9 shows that at 

zero crossflow velocity and 100 percent of design speed, the calcu- 

2 

lated value of total thrust loading is 809 lb/ft as compared to 

2 

the observed value of 819 lb/ft . At 70 percent of design speed 
. (Fig. 10) the calculated value of total thrust loading is 405 
lb/ft as compared to the observed value of 411 lb/ft . The data 
also shows that the calculated values show the same trends with in- 
creasing crossflow velocity as the observed values. 

The data on Figures ll and 12 show that at zero crossflow ve- 
locity the calculated and observed values of axial-force loading are 
almost identical. Over the 1 crossflow range tested, there is good 
agreement between the calculated and observed values, except at 
zero angle of attack and 100 percent of design speed test condi- 
tions. At these test conditions, the observed values appear to 



Total thrust loading, F/6A, lb/ft 



Total thrust loading, F/6A, lb/ft 



Figure 10. - Comparison of Computed Values from Correlation Equ- 
ation. and Observed- Values of Total Thrust Loading at 70 Percent 
of Design Speed. 



Axial -force loading, F /6A, lb/ft 


«1' 



Crossflow velocity, V, ft/sec 

Figure 12. - Comparison of Computed Values from Correlation 
Equation and Observed Values of Axial-Force Loading at 70 
Percent of Design Speed. . 
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have more scatter than at the other test conditions. However, no 

reason could be found for this apparent deviation. 

The values calculated from the correlation equation will now 

be used to discuss the variation in loadipg parameters with the 

various test variables. The data on Figures 9 and 11 show that at 

zero crossflow velocity for 100 percent of design speed, the total 

2 

thrust loading is 809 lb/ft .as compared to the axial-force load- 

. . . 2 • ■ ' . ■ 

ing value of 744 lb /ft . At 70 percent of design, speed (Figs. 10 

• ■ • 2 
.and 12), the total thrust loading is 405 lb/ft as compared to the 

2 

axial-force loading of 350 lb/ft . • This observed difference can . 
be attributed to the fact that a portion of the fan momentum thrust 
acts on the wing because the inlet bellmouth of the lift-fan simu- 
lator is not of an infinite diameter. If the inlet had an infinite 
diameter, then the maximum lift force generated would be equal to 
the exit momentum thrust less the fan base force. Consequently, in 
this experiment, thp full momentum thrust will not be. measured on 
the load cell force balance resulting in a reduced value of axial - 
force loading. 

A comparison of the data on Figures .9 thru 12 show that the 
effect of increasing crossflow velocity is not the same for total 
thrust loading as for axial force loading. At 100 percent of de- 
sign speed and zero angle of attack (Fig. 9), the total thrust 

2 

loading decreases uniformly from 809 lb/ft at zero crossflow 

2 

velocity to about 655 lb/ft at the highest crossflow velocity 

of 240 ft/sec. However, the axial-force loading (Fig. 11) 

2 

tends, to decrease slightly (from. 744 Ib/ft at zero crossflow 



39 - 


velocity) with increasing crossflow velocity to minimum value of 
2 

696 lb/ft at about 150 ft/sec. Further increases in crossflow 
velocity to 240 ft/sec cause an increase in axial-force loading ap- 
proaching the same level as that obtained at zero crossflow 
velocity. 

At 70 percent of design speed (Fig. 10), the total thrust 

loading decreased slightly with increasing velocity. However, the 

axial-force loading (Fig. 12) increased uniformly with increasing 

crossflow velocity. For example, at zero angle of attack, the 

• 2 

axial-force loading increased from 350 lb/ft at zero crossflow 

2 

velocity to a value of 420 lb/ft at the highest crossflow veloc- 
ity of 240 ft/sec. 

‘ The observed difference in trends between total thrust load- 
ing and axial-force loading with increasing crossflow velocity may 
be partially explained by the fact that crossflow velocity may ef- 
fect the total- thrust loading .in two ways.' First, crossflow ve- 
locity tends to cause flow separation in the fan inlet. Inlet . 
flow separation was actually observed in this experiment. When 
flow- separation occurs it reduces the effective flow area of the 
fan which reduces the weight flow through the fan. The inlet area 
blocked by flow separation increased with increasing crossflow ve- 
locity and increased angles of attack. Inlet flow separation was 
also found to be a function of corrected tip speed since more in- 
let flow separation was observed at the design tip speed (980 
ft/sec) than at 70 percent of design tip speed. 
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The second effect of crossflow is to produce flow distortions 
at the inlet and within the fan flow passage. This effect reduces 
the overall total pressure at the exit of the fan. Both of these 
effects reduce the total thrust loading, and can explain the rea- 
son for the observed decrease in total thrust loading with in- 
creasing crossflow velocity and increasing angle of attack. 

' A comparison of the correlation equation and the observed 

values of the force loading difference is shown in Figures 13 and 

14. The data show that at zero crossflow velocity the total 

thrust loading is greater than the axial force loading at both 100 

and 70 percent of design speed. The correlation equation shows 

that at 100 percent of design speed (Fig. 13), the difference is 

2 

-65 lb/ft ; and at 70 percent of design speed (Fig. 14) the dif- 

2 . 

ference is -58 lb/ft . The reason for the difference has been ex- 
plained by- the previous discussion. • 

The data in Figures 13 and 14 show that the force loading 
difference increases uniformly with increasing crossflow velocity, 
and that the force loading difference increases more rapidly at 
.+10° than at -10^ angle of attack. 

A comparison of the correlation equation and observed data 
for total thrust loading as a function of axial-force loading is 
presented in Figure 15. The data showed that at zero crossflow 
velocity, the observed data is displaced by about a constant value 
from the correlation equation. The effect of angle of attack aind 
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Angle 



Crossflow velocity, V, ft/sec 

Figure 13. - Comparison of Computed Values from Correlation 
Equation and Observed Values of Force Loading Difference 
at 100 Percent of Design Speed. 



Force loading difference AF/6A, lb/ft- 
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Figure 14. - Comparison. Of Computed Values from Correlation 
Equation and Observed Values of Force Loading' Difference 
at 70 Percent of Design Speed. 



Total thrust loading, F/6A, lb/ft 
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crossflow velocity is to vary the observed data from the linear 
correlation equation shown. Further examination of the data on 
Figure 15 indicates that a family of lines with the same slope as 
the original correlation equation could be drawn through the data 
to represent constant values of angle of attack or crossflow ve- 
locity. The displacement of the lines from the original correla- 
tion equation would vary with . crossflow velocity and angle of at- 
tack. ' It may then be possible to determine correlation equations 
for the family of lines to define the performance characteristics 
pf the lift-fan simulator for the full range of dependent vari- 


ables . 



6.0 CONCLUSIONS AND AREAS FOR FURTHER RESEARCH 

From the results of this experiment, the following major con- 
clusions were reached: 

The multiple linear regression analysis technique may be used 
to obtain significant correlation equations for the dependent vari- 
ables: total thrust loading, axial-force loading, and the force 

loading difference .as a function of the independent test variables : 
crossflow velocity, corrected tip speed, and wing angle of attack. 
The results of the multiple linear, regression analysis showed that 
for all cases the correlation coefficient exceeded 0.988 which in- 
dicates a near perfect correlation. The F values obtained from 
the F-tesf for testing the null hypothesis '’Hq':' that' .there is no 
correlation, were greater than the critical value of the F- 
distributions for a 0.01 probability level. Therefore, the null 
hypothesis H^: that there is no correlation., was rejected and the 

alternative hypothesis, H^: that there is ;a significant correlay 

tion, was accepted. 

The results of the linear regression analysis also showed that 
there is a significant correlation for the total thrust loading as 
a function of axial-force loading. 

Since the linear regression analysis technique can calculate 
reliable prediction model equations, it would be- possible to reduce 
the amount of instrumentation required to. measure lift-fan simula- 
tor performance. For example, only a load cell force balance and 
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limited pressure instrumentation would be used to measure the lift- 
fan axial force, and the prediction model equation would then be 
used to calculate the total thrust load. Or if the lift-fan simu- 
lator design could not accommodate a force balance, the prediction 
model equation could be used to calculate the axial-force loading. 

The prediction model equation was also useful in determining 
the relationship between the axial-force loading and the calculated 
’total thr.ust loading. Crossflow velocity was found to be. the most 
significant variable in the force-loading difference correlation 
equation. The correlation equations were also used to smooth the- 
observed data. Smoothing the data helps to visualize the interfer- 
ence effects of the fan flow and external airflow (crossflow veloc- 
ity) on the lift-fan simulator performance. It also helps to ex- 
plain the reasons for the differences between the axial-force load- 
ing and the total thrust loading. 

The experimental test results showed that the total thrust 
loading decreased uniformly with increasing crossflow velocity be- 
cause of inlet flow separation and flow distortion effects on the 
exit total pressure. The axia^force loading tended to first de- 
crease slightly and then increase with increasing crossflow- veloc- 
ity. This effect is a result of the combination of decreasing 
total thrust loading and increasing circulation pressure forces 
as the crossflow velocity increases. 

The results of this work indicate that prediction modeling 
techniques could be valuable in designing and analyzing lift-fan 
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simulator tests. It could be used to reduce the amount of instru- 
mentation required and to reduce the number of test runs required 
to define the performance characteristics of lift-fan simulators. 

The use of this technique could reduce the cost of building and 
testing lift-fan simulators. 

However, a large amount of work remains to be done before the 
method presented in. this thesis c'puld actually be used routinely 
for any VTOL -lift-fan test programs, The technique would have to 
be applied to other lift-fan simulator test programs to confirm the 
findings of this work. As previously stated, three other lift-fan 
simulator tests have recently been completed. This technique may 
also be applied to the results of those lift-fah simulator tests. 

In addition, the techniques presented here should be applied 
to the test results of the lift-fan simulator tests with the' exit 
louvers installed at the exit of the' fan. The effect of exit lou- 
vers on the .fan performance characteristics of. the lift-fan simula- 
tor is a fourth independent variable- to be included in the correla- 
tion equation. 

The. goal of this thesis was to develop a prediction technique 
to help in determining the performance characteristics of lift-fan 
simulators. The technique has tremendous potential, and it is hoped 
that this work will be a useful tool in the design and analysis of 
future lift-fan simulator test programs. 
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APPENDIX A 


METHOD OF CALCULATING TOTAL THRUST 


This chapter describes the methods used in calculating the fan 
and turbine momentum thrust, axial force, and fan base force from 
the measured state point values of pressure, temperature and air- 
flow angles, and axial force recorded at the. various test points. 

‘ A. 1 Momentum Thrust : . ' • 

The fan and turbine momentum thrust are . calculated by the 
same procedure. The momentum thrust is calculated from the meas-. 
ured weight flow, and the ideal exit axial velocity calculated from 
the ideal expansion from a calculated total pressure to an undis^- 
turbed ambient pressure in conjunction with a thrust coefficient 
for outlet duct losses. 


Fan Momentum Thrust . - The equation for the fan momentum thrust is 



w V C • 
1 ze tf 


lb 


(Al) 


where C^ is the thrust coefficient determined from experimental 

data of similar exit ducts. For this application • 1.0. The 

2 

acceleration due to gravity, g = 32.2 ft/sec . V is the ideal 

Z6 

exit axial velocity and is deterniined from the following equation. 
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where p Q is the undisturbed ambient static pressure measured in 
the test section, P is the average total pressure at the fan exit 

© v 

weighted according to the distribution of mass flow in the fan pas- 
sage (such terms are denoted as mass averaged terms); is the 

average airflow angle; is the mass averaged total temperature, 

and y is the ratio of specific heat of air,, y = 1.4. 

•The weight flow w^ is calculated from-. the correlation 

Wj^/e/o = lb/ sec (A3) 

where p^ is the' fan inlet, static, pressure 3/4 of an inch above 
the rotor leading edge. The correlation is obtained from a compu- 
ter program based on the potential flow theory of referenpe 22 for 
an inlet sufficiently far upstream pf the rotor to be considered 
unaffected by the rotor. The weight flow is also adjusted for the 
effects of flpw separation on the inlet bellmouth. 


Turbine -Momentum Thrust . r-.The equation for turbine momentum thrust 

t 

is . 


F = 
T 


W T V zT C tf 

g 


lb 


(A4) 


where the thrust coefficient C = 0.98. V _ is the ideal exit 

tf zT 

axial Velocity and is determined fyom the following equation. 




(A5) 


1 


, ft/sec 
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Where is the mass averaged total pressure at the turbine $xit, 

IN 

T^ is the mass averaged total temperature, and 3^ is the average 
airflow angle. The weight flow w^, is measured with a calibrated 
flat-plate orifice. 

A, 2 Fan Base Force 

. The fan base force Fg is calculated from the measured- static 
pressure p^ acting On the lower. surface of the fan hub. 

.. . f b - <p„ - w lb - (A6) 

Where p is the ambient static pressure in the test section, and 

2 

Ag is the area on which p^ acts (A^ = 0.227 ft for this exper- 
iment). 

A. 3 Axial Force 

The axial force is measured on the load cell force balance. 

The force balance consists of three 500 lb load- ceils equally 
spaced on a 22.5-in. diameter circle (Fig. 3). The total axial 
force is the sum of the three load cell readings times a calibra- 
tion constant . 

F lc = (F 1 '+ F 2 + F 3 )K f ,.lb. (A7) 

where F^ , F 2 , and F^ are the individual load cell readings . 
is a calibration constant determined from calibration tests of the 
load cell force balance. 



APPENDIX B 


DESCRIPTION OF MULTIPLE LINEAR REGRESSION COMPUTER PROGRAM 
27 

RAPIER is a FORTRAN IV program for multiple linear regression 
analysis. RAPIER computes the variance-convariance matrix of the 
independent variables, regression coefficients, t-statistics for in- 
dividual tests, and analysis of variance tables for overall testing 
of regression. There is a choice of three strategies for the vari- 
ance estimate to be used in computing t-statistics. A backward re- 
jection option method, based on. the first dependent variable, may be 
used. A critical significance level is supplied as input. The 
least significant independent variable is deleted and the regression 
recomputed.. This process is repeated until all remaining variables 
have significantly nonzero coefficients. 

The model equation is assumed to be of the form 

y . = b 1 x. + b„x. +'••••+ b T x. _ + e. i = 1, . . .N (Bl) 

. ■'i 1 i^ 2 ± 2 . J iJ i 

where: 

y^ is the response variable. 

x is the transformation of the independent variables. 

1J 

The. subscript i indicates the values associated with the 

i^ observation. 

The subscript J denotes the number of Independent vari- 
ables observed. 

■N .is the number of observations. 

e^ denotes the difference between the observed value and 

the expected value of y^. 
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The program then computes the regression coefficients b^, b^,.... 

b^; the (Analysis of Variance) ANOVA table, the coefficient of de- 
2 

termination r , and the coefficient of correlation r; the stand- 
ard error of the regression estimate. In addition an F test is 
made for testing the null hypothesis H^: that there is no correla- 
tion, that all the correlation coefficients are zero, the alterna- ' 
tive hypothesis is ..H^: that there is a significant correlation . 

The F test is 

p - Mean squares of the regression 
Mean squares .of the residual 

If the calculated value of* F is greater than the critical value of 
the F distributions for the number of degrees of freedom, then the 
hypothesis Hq is invalid indicating that there is a significant 
correlation. and the regression coefficients are not zero. 

If the backward rejection option of the program is used a 
t-test .is performed on each regression coefficient. The term cor- 
responding to the minimum t is dropped and the regression is re- 
computed. This process is repeated until all remaining coefficients 
are significant at the specified level of confidence. 

Then for. each observation (input data) the observed response 
(y) , the calculated response (y calc) ; the residual (Y obs - Y calc = 
Diff) , and the studeqtized residual z = ^ is computed and 

printed out. Two other statistics which may be used to test the 
normality of the empirical distribution, skewness and kurtosis, are 
calculated. Skewness should be nearly zero and kurtosis should.be 


nearly 3. 
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The program has a prediction option, in which the regression 
equations are used to predict an estimated response at some condi- 
tion of the independent variables. The program accepts input vec- 
tors x and computes the predicted response, y; the variance of the 
regression line, the standard deviation of regression, variance of 
the predicted value and the standard deviation of the predicted 
value . ... . • ; 

- The • program also has a computer plot routine which plots the 
residuals for the dependent variable versus. the independent vari- • 
ables and the .residuals for the dependent variable versus the pre- 


dicted value. 
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*b 

H- 

C tf 

F.'- 

F B 


F. • 
x 


LC 

F T 

F/6A 

V 6A 

AF/6A 

H 0 

H 1 

P 

P A . 


GLOSSARY OF SYMBOLS 

*> 

2 

fan frontal area, 1.263 ft 
fan hub base area, 0.227 ft^ . 
regression coefficient 
thrust coefficient 

total thrust force, lb . 

fan base force', lb 
fan momentum thrust, lb 
individual load Cell reading, lb 

axial-force measured on load cell force balance, lb 
turbine momentum thrust, lb 

2 

total thrust loading parameter, lb/ft 

2 

axial-force loading parameter, lb/ft 

2 

force loading difference, lb/ft 
statistical hypothesis to be tested 


s • 

P 1 

AP 


alternative hypothesis to be accepted if Hq is judged 
to be false 

2 

total pressure measurement, lb/ft 

• 2 

total pressure measurement for airflow angle, lb/ft 

2 

mass averaged fan exit total pressure, lb/ft 

. •• 2 

static pressure on lower surface of fan hub, lb/ft 

2 

static pressure measurement, lb/ft 

2 

fan inlet static pressure, lb/ft 

2 

differential pressure measurement, lb/ft 
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T total temperature measurement, °R 

A Q_ 

T g mass averaged fan exit total temperature, ^ 

T^, mass averaged turbine exit total temperature, °R 

AT differential total temperature measurement, °R 

U T //a fan corrected tip speed, ft/sec 

V crossflow velocity (free stream airflow) , ft/sec 

V ideal exit axial velocity of fan, ft/sec 

V ^ ideal exit axial velocity of turbine, ft/sec 

w^ fan weight flow, lb/sec 

w^ turbine- weight flow, Ib/sec 

% 

y dependent variable for regression analysis . 

.a wing angle of attack., deg 

3 airflow angle, deg 

6 the ratio of inlet total pressure to standard sea-level 

.pressure of 2116.2 lb/ft^ 

. 0 the ratio, of inlet total pressure to standard sea-level 

- temperature 518.7° R . 

y • ratio of the specific heat of air, y = 1.4 . 



